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Abstract 

We study the spectrum of an integrable antiferromagnetic Hamiltonian of the gl(M|A'^) 
spin chain of alternating fundamental and dual representations. After extensive numer- 
ical analysis, we identify the vacuum and low lying excitations and with this knowledge 
perform the continuum limit, while keeping a finite gap. All gl(n + A'^IA'^) spin chains 
with n,N > are shown to possess in the continuum limit 2n — 2 multiplets of massive 
particles which scatter with gl(n) Gross-Neveu like S'-matrices, namely their eigenvalues 
do not depend on A*'. We argue that the continuum theory is the gl(Af|A) Gross-Neveu 
model. We then look for remaining particles in the gl(2m|l) chains. The results suggest 
there is a continuum of such particles, which in order to be fully understood require 
finite volume calculations. 
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1 Introduction 

At present, the continuum limit of antiferromagnetic Lie superalgebra spin chains is poorly 
understood. High degeneracies or continuous spectra |EFS05[I5S07[|IJS08] and singular IR 
scattering amplitude behavior |EFS05IISP10) emerge in the continuum large volume limit. We 
interpret this as a hint of new physical phenomena. Strikingly, when g-deformed, connections 
with long known puzzling theories with no obvious supergroup symmetry appear [SPlOj . 

Previous, above mentioned, research has mostly dealt with osp(2|2) ~ sl(2|l) spin chains. 
We adopt a general approach, taking as a model |EFS05| . and hoping that the greater 
gl(M|A^) symmetry and the relatively new cohomology techniques developed in [CCMSIO] 
will allow us to better understand the nature of these new phenomena, thus, suggesting how 
to correctly interpret, treat and, eventually, apply them. The goal is to ultimately formulate 
a meaningful factorizable scattering theory for the continuum limit of superspin chains. This 
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is an even harder task then it sounds, because there are few massive relativistic field theories 
with supergroup symmetry that have been understood, at least partially |SWK02llSP10l 
iBLnOllOLLOO] . 

So, clarifying the continuum limit of superspin chains must ultimately result in a better 
understanding of formally simple field theories such as the gl(M|iV) and osp(i?|2S') Gross- 
Neveu models, which play an important role in disordered electronic system [BerilGLLOOl 



ILeCOQ) . provide instances of continuous families of conformal field theories |CS09bl CMQ+10 



and appear in strong- weak coupling dualities with sigma models |CS09bl[GLL00[[C"R09) . Due 
to strong violation of unitarity, it is not obvious how to treat these super Gross-Neveu models 
directly in the continuum by standard bootstrap methods |ZZ79[[ORW87] . Insight from 
coordinate Bethe ansatz for chiral Gross-Neveu models |AL80allAL80bllAL79[|DdV89j has led 
us to consider the spin chain of gl(M|iV) alternating fundamental and dual representations 
as the most natural candidate for an integrable discretization of the gl(M|A^) Gross-Neveu 
model. The spin chain of gl(M| A^) fundamental representations is discarded because it cannot 
lead to a relativistic field theory |Sal001lEK94| . This can be immediately seen from the fact 
that the dual of any g\{M\N) representation appearing in the tensor product of gl(M|A^) 
fundamental representations does not belong to this tensor product if A > 0. 

The article is organized as follows. In sec. [5] we define the spin chain and its integrable 
dynamics. In sec. |3] we explain how to perform numerical calculations for spin chains of 
modest size, but arbitrary gl(M |A^) symmetry, using the walled Brauer algebra. Then we 
present numerical results on the spectrum. Sec. |4] is of primary technical importance. We 
discover a class of solutions to the g\{M\N) Bethe ansatz equations (BAE) which can be 
fully characterized in terms of solutions to the gl(A/ — 1\N — 1) BAE. This provides an 
explicit embedding of the gl(M — 1|A^ — 1) spectrum into the gl(M|A^) spectrum. We give 
an algebraic explanation to this relationship. Further, in sec. [S] we identify the vacuum 
solution of gl(M|A^) BAE and write it down explicitly in terms of the gl(M — N) vacuum 
solution. We then classify gl(M|A^) excitations that lead to a gl(M — A^) like spectrum and 
the simplest ones that do not. For the latter, we analyze numerically the form of the solutions 
to the BAE. Finally, in sec. [5] we consider the continuum limit of spin chains with alternating 
inhomogeneities. Excitations which can be characterized in terms of solutions to gl(A/ — N) 
BAE lead to a spectrum of massive particles with gl(M — N) Gross-Neveu mass ratios and 
S'-matrix eigenvalues. We argue that the continuum limit of the spin chain is the gl(M|A'') 
Gross-Neveu model. We end with analyzing some low lying excitations of the gl(2TO|l) spin 
chain which are not of gl(2m — 1) type and lead to new massive particles. 

2 Transfer matrices, Hamiltonians and their spectra 

Let V be the fundamental module of gl(M|A), V* denote its dual and p : gl(M|A^) ^ gl{V) 
together with p : gl(M|A^) ^ g^iV*) be the corresponding representations. If {e^}*^^ 
is a graded basis of V with grading |a| :— \ea\ G Z/2Z and {e^j^^Ji^ is the dual basis 
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e"(e^) — then Eafi are the standard generators of gl(M|A^) acting m V as Ea^ ■ = 
8fi^ea and in V* as E^fi ■ e' = -(-I)(l"l+I''l'l"l5a7e'^. Let us label the V factors of the 
spin chain C{L) = {V ® y*'^®^ fj-Qj^ i^f^ j-igj^t by a subscript 1, . . . , L. Similarly, we label 
the V* factors by a subscript 1,...,L. For £ e EndF we denote by £k G EndC(L) the 
endomorphisms acting as £ on Vk and trivially, up to grading signs, everywhere else in the 
chain. Similarly, for £ G EndF*, fj. £ EndC(i) will act as £ on and trivially, up to 
grading signs, everywhere else. 

The C(L)-endomorphisms Pki = (—1)^^^ Pk{Eap)pi{Ei3a) provide a representation for the 
symmetric group acting on the V factors of C{L). Similarly, Pj^j = (—1)^^^ p^.{Eap)pi{Ei3a) 
generate a representation for the symmetric group acting on the V* factors of C(L). On 
the other hand, Q/^i = —{—iy^^pk{Eap)pi{Epa) generate a representation of the Temperley- 
Lieb algebra T2L{n) with loop weight n = M — N. Together, the P's and Q's generate the 
gl(M|A^)-centralizer algebra of C{L), that is the set of all endomorphisms of the spin chain 
that commute with the gl(M|A^) action [SerOlj . This centralizer algebra is a representation 
of the walled Brauer algebra BL,L{'n), which can be viewed as a subalgebra of the Brauer 
algebra B2L{n). We shall discuss in detail the algebra Bl^lIti) and its representations in 
sec. 13.11 

In terms of walled Brauer algebra generators we introduce the i?-matrices 

Rij {u) + Pij R^j (u) = u - Q^j 

Rij{u) = u - Qij ^ u + Pij , 

which satisfy the Yang-Baxter algebra 

Rij{u - v)Rik{u)Rjk{v) = Rjk{v)Rik{u)Rij{u - v) (2.1) 
R^jiu - v)R,k{u)Rjk{v) = Rjkiv)Rikiu)Rtjiu - v) 
Rij{u -v + n)Rik{u)Rjk{v) = Rjk{v)Rik{u)Rij{u - v + n) 

and similar "dual" relations, that is the above relations with all barred indices replaced 
with unbarred ones and all unbarred indices with barred ones. We define two one-parameter 
families of monodromies 

Ta{u) = R,i{u + n/2)RaL{u) . . . R,i{u + n/2)Rai{u) (2.2) 
fa{u) = R-^i{u)RaL{u + n/2) . . . R^i{u)Rai{u + n/2) (2.3) 

acting on Va ® C{L) and Va <8> C(i), respectively, and corresponding transfer matrices 

t{u) = straT,(u), t{u) = straTa(u) . (2.4) 

Yang-Baxter relations (j2.1l) imply a Yangian structure given by the two relations 

Rabiu - v)Ta{u)Tb{v) = Tb{v)Ta{u)Rab{u - v) (2.5) 

i?„5(u - V + nl2)Ta{u)fb{v) = fb{v)Ta{u)R,i{u - V + n/2) (2.6) 
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and their duals. The fohowing commutation relations immediately follow 

[t{u),t{v)] = [t{u),i{v)] = [t{u),i{v)] = . 

The nested algebraic Bethe ansatz for the most general gl(M|7V) spin chain was considered 
in IBROSj . The Bethe ansatz equations and the spectrum of t{u) formally depend on the 
nesting order, that is an ordering of a basis of V and the induced ordering of the dual basis 
of V* . If the basis {ea}a^Ji^ of V diagonalizes the Cartan subalgebra, then, without loss of 
generality, we label the basis vectors so that the total ordering reads 

ei > 62 > • • • > CM+N , (2.7) 

where, however, we keep unspecified the grading of the basis vectors. Let wt{ea) — ea 
denote the weights of basis vectors called fundamental weights. The ordering (12.71) induces 
a weight space ordering ei > • ■ • > cm+n which fixes the simple root system Aq = {a^ := 
Ci — ei+i}fl^^~"^. The choice of grading, which we denote by S = {ctq = \ea\ = {aW^^Ji'^ , 
is equivalent to the choice of a Cartan matrix, or a Dynkin diagram. Changing the grading 
can be equivalently seen as changing the total ordering (|2.7p . As a result, the simple root 
system, the positive root system and the Borel subalgebra changes with E. So, keep in mind 
that the notion of highest weight always depends on the grading choice and changes when E 
changes. 

Define the operator matrix elements of the monodromies (12.21 12. 3p as T = Eij ® Tij , where 
Tij G EndC(L) and T — Ta, Tg. Choosing the reference state to be the highest weight state 
of C{L), the eigenvalues of t{u) can be written in terms of polynomials {Ta)ii{u)VL = Ki{u)VL 

r(w + (-l)|i|)^(?i + n/2)^, 1 = 1 
Aj(u) = < u^(M + n/2)^, 2<i<r 

[w^(u + n/2- (-1)1^-^+^1)^, i = M + N 

and simple root Q-polynomials 

Qk{u) = ll{u-uf) 

as follows 

^ Qk-i[u)Qk[u) 

(k) 

where the Uj , k = 1, . . . ,r = M + N — I are the Bethe roots appearing at the fc-th step of 
the nesting. There are solutions of the following system of nested Bethe ansatz equations 

Ak+i{uf^) Q,_i(uf + (-l)W)Q,(uf -(-l)l'=l)g,+i(uf ) 

ensuring the analyticity of eigenvalues (j2.8p . where k = 1, . . . , r and j = 1, . . . , v^''"^ and it has 
to be understood that Qq{u) = Qm+n{u) = 1. We stress that the BAE (|2.9p are equivalent 
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to the analyticity requirement of the transfer matrix (j2.8p if and only if none of the Bethe 
roots of the same type coincide, which is an essential requirement in the algebraic Bethe 
ansatz construction. 

The BAE take a more familiar look }OW86) 

L v" 

ne(A.,a)(4"^-j/fc)e(A,..)(4"'-y^)--(-i)i"i n -A^^) ' (2-10) 

k=l ^6Ao i=l 

when written in terms of new variables 



X 



3 2 ^ 

l=\ 



the Takahashi functions 



et(a;) 



(-l)l'l , (2.11) 
i</2 



- it/2 ' 

the highest weights A^^; = ei and Aj — —^m+n of modules Vk and in C(L), and the weight 
space scalar product {ei,ej) — 1)'*'. The degree of a root ak = — f-k+i is defined 
as jafcl = + j/s + 1|. We shall be mostly considering the homogeneous case (|2.2[ 12. 3p 
corresponding to = j/j, = 0, although the inhomogeneous deformation yk,yk ^ shall also 
be required. There are j = 1, . . . , z^" equations for every a £ Aq. 

The weight of the reference state is wt(f2) — J2k=i ^k + J2jc=i ^k = -^(^i ^ ^m+n)- Bethe 
vectors lo described by the Bethe roots (|2.10p are highest weight vectors of weight 



du 



wt(w) = wt(^2) - ^ afcz/^'^) . (2.12) 
fc=i 

We define the dynamics of the spin chain by the momentum and Hamiltonian operators 

_o °SAi(u)Am+^(^) ' ^^P*^-( 1) Ai(0)Am+^(0) • ^'-'-^^ 

In order to have explicit expression for the spectrum of these operators, the eigenvalues (|2.8p 
of t(u) are not enough. One also needs to evaluate the eigenvalue of i{u) on a given Bethe 
eigenvector of t{u). A fundamental difference w.r.t. gl(A^) spin chains is that one cannot solve 
this problem by fusion. This is because tensor products of the fundamental representation 
V of sl(Af |A^) will never generate the dual representation V* as a direct summand, nor even 
as a subquotient. To develop a clear idea about how this should be done, let us recall that 
a set of Bethe vectors for t{u) can be constructed in the framework of the algebraic Bethe 
ansatz (ABA) by using only the commutation relations (|2.5p . Then, the eigenvalue of t{u) on 
such a Bethe vector can be calculated, in principle, by using the second type of commutation 
relations (|2.6p . We shall not pursue this rather tedious route. Instead, we guess the eigenvalue 
of t{u) on a given Bethe vector of t{u) as follows. 

First, notice that a different set of Bethe vectors can be obtained by performing the 
ABA for i(w), that is by using the commutation relations dual to eq. (j2.5p . We perform the 



5 



nesting by ordering the dual basis vectors {e"}*^^ of the auxiliary space V* according to 
their weights wt(e") — — £„ 

e^i+N > . . . > e2 > gi _ (2.14) 

We denote the Bethe roots appearing at step k of the nested ABA by u''j^'^~^^~''\ because the 
simple root which must be associated to them is wt(e*^"'"^~'^"'"^) — wt(e*^+^^'^) — au+N-k- 
Then, the eigenvalues of t{u) can be written in term of polynomials {Ta)ii{u)Q. — Ki(u)U, 

{(M-(-l)|i|+n/2)^M^, i^l 
{u + n/2)^u^, 2<i<r 
(M + n/2)^(u+ (-1)1^^+^!)^, i^M + N 



(fe) 



and simple root polynomials Qk{u) — Y\'^j^i{u — Uj ) in complete analogy with 

A(„).'x:Vi)'-'A.w '^'--'''-'-^'''''g';v'-''''r (2.5, 

^ Qk-i(u)Qk(u) 



k=l 



Analyticity conditions for A(u) take the same form as eqs. (|2.10p in terms of variables 



. M+N 

{ E (-1)'" (2-16) 

l=k+l 

and parameters A Bethe vector constructed in this way, which we denote by w, has 

weight 

r 

Wt(w) = Wt(0) — UkV^^^ . 

k=l 

At this point, the eigenvalues of t{u) w.r.t. the second set of Bethe vectors ui is not 
known. If we can match Bethe vectors uj with Bethe vectors Co then the problem is solved. 
Due to the subtle completeness issue of ABA for super spin chains, it is not clear at all if 
the matching can actually be performed. We assume it can and we do it as follows: a Bethe 
vector Lu of t{u) is identified with a Bethe vector ui of t{u) if 

^ 1 rJfc) _r^WiP<'=' fc = l---^- (2-17) 

The first condition in the braces means wt(tj) — wt(a)). Eqs. (|2.8[ 12.151 12.111 [2.161) allow us 
to compute both the eigenvalues of t{u) and t{u) on a Bethe vector a; = w in (|2.17p . 

After this long detour, we come back to the spectrum of the operators (j2.13p . which can 
now be explicitly computed 

E = -y24^ E^T-r (2-18) 

P = E(-l)'"^i(4'') + E(-l)"'^'^'^i(4''^) i^od 27T , (2.19) 
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where 9t{u) — iloget{u) + tt = 2tan~-'^ ^ with some fixed branch. First thing to be noticed 
is the exphcit dependence of the definitions (j2.13p on the grading. Therefore, it iooks like 
the type of the chain - ferromagnetic or antiferromagnetic - also depends on it, which is 
also suggested by the grading signs in (|2.18l) . However, the charges H and P can be written 
explicitly as a representation of an element of the periodic walled Brauer algebra 

H = f2 -(-!)'" - (-1)!'"+^' + +P—-1 {{Qu+i^Q^^ + {Q—u,Qu}) (2.20) 

where all (un)barred indices are defined modulo L and the afhne generators are expressed 
in terms of non-periodic walled Brauer algebra elements Pli ■= Pil ~ P12 ■ ■ ■ Pl-il ■ ■ ■ P12, 
Pjj ■■= Pji = • • ■ Pi~=Tl-- ■ Pvl and Qii := Qii = PilQllPil- The following general 
formulas have been used 



dlogt(O) A 



«5'(f)«..(i)+E«;'(5)A«(o)*(i) 



du ^ 

i=l i=l 
i—l i—l 

The latter can be derived using only the cyclicity of the supertrace, which holds for a graded 
tensor product when even endomorphisms are considered, relations of the form PijRjk{u) = 
Rik{u)Pij, PijRjj^{u) = R^^{u)Pij, stTaPaj = 1 and their duals. We see from eq. (|2.20p 
that the grading enters in the definition of the Hamiltonian only as a shift, therefore having 
nothing to do with ferro or antiferromagnetism. The equivalence of the solutions of BAE 
in different gradings, and therefore of spec 77, is at present somewhat understood in terms 
of particle hole transformations |Tsu98| , although the relationship between Bethe vectors in 
different gradings not at all. 

A relation between the spectrum of H and —H can be constructed, but one has to consider 
different chains. Let Hm\n denote the integrable Hamiltonian of the gl(M|iV) spin chain. 
Then on has the following relation 

Hm\n = ^Hn\m ■ (2-21) 

To prove it one uses the algebra homomorphism between the walled Brauer algebras BL j^^n) 
and BL^L{—n) provided by P M> — P and Q h- !■ —Q. This automorphism is realized in the spin 
chain representation by the shift of the grading function |z| i-> |i| + 1, which maps gl(M| A^) n- 
gl(A'^IM). As we shall see in the next section, the sign in front of our Hamiltonian ((S^ 
ensures that we are dealing with antiferromagnetic spin chains for n = M — N > 0. Eq. (|2.2ip 
allows to fold back the gl(A^|Af) spin chains with A^ < Af to gl(Ajf|A^) spin chains with 
n > by changing the sign of the Hamiltonian (|2.20p . but now they will be ferromagnetic. 
As the gl(A|A^) chain in eq. (|2.20p is poorly defined, we shall discard it from the main 
discussion and come back to it only in the conclusions. From now on we restrict to gl(A/|A^) 
antiferromagnetic spin chains (|2.20p for which n = M — N > 0. 
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I I LIKLI LIXLI LIxLI 

I ... L I ■■■ i ... L I ■■■ i ... L 1 ••• i ••• L 

I Ei Pi Pi_iPjPj+i 

Figure 1: The identity I and the generators E^, P^ of the Brauer algebra of dimension (2L — 1)!! 
are represented on the left; the walled Brauer algebra generator P^-iPPi+i is represented on 
the right. 

3 Numerical diagonalization of the Hamiltonian 

Diagonalizing the matrix (I2.20p is not the smartest thing to do if one is solely interested in 
its spectrum, especially if one intends to consider on the same footing all the gl(M| A^)-spin 
chains. This is because many eigenvalues have multiplicities corresponding to the dimension 
of irreducible gl(MjA^) representations appearing (generically as quotients of submodules) in 
the spin chain. These multiplicities quickly grow with the rank and the computing power 
per eigenvalue increases respectively. An approach which allows to select the eigenvalues 
corresponding to a given gl(M|iV)-irreducible representation and eliminate the corresponding 
degeneracy would be considerably more efficient. To implement this approach one interprets 
the Hamiltonian H in eq. (I2.20p as an element H of some algebra, namely the walled Brauer 
algebra B^xin — M—N), in a particular representation provided by the gl(M|A^)-centralizer 
of the C{L) spin chain. As we shall explain shortly, the algebra BLx{n) can be abstractly 
defined independently of the gl(n + N\N) spin chains and, in particular, it does not depend 
on N . Ccntralizers of gl(n + N\N) spin chains provide N dependent representations for 
BL,L{n). The important thing is that the representation theory of BLx{n) is understood 
well enough, so that one can find the spectrum of the algebraic Hamiltonian H by working 
in other representations where its spectrum is much less degenerate compared to that in the 
gl(n + N\N) spin chain, namely it does not depend on TV. 

3.1 Walled Brauer algebra and its standard modules 

The walled Brauer algebra BL^in) can be conveniently viewed as a subalgebra of the Brauer 
algebra i?2L("-)- The defining relations of B2L{n) can be found in jGCA98| . The words of 
B2L {n) admit a representation as graphs on AL labeled vertices with 2L edges connecting the 
vertices pairwise in all (4i — 1)!! possible ways (crossings are allowed). The identity I of the 
Brauer algebra and the generators E^, P^ are represented by the graphs on the left in fig.[Il 
In order to multiply the diagrams one arranges the first 2L vertices horizontally with 
the remaining 2L vertices on top of the first ones. The product of a diagram di with a 
diagram 6,2 is the diagram did2 obtained by i) placing the diagram di on top of the diagram 
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• •(giv 




• ^^•®((2,l),(l,2)).v 



Figure 2: Example of diagrammatic multiplication of a basis element of A3 3(A, /x), where 
X,fj,\~ 2 and v is an element of the module S{X) x S{fi) of Sym(2) x Sym(2). 

d2, ii) identifying the top of the diagram d2 with the bottom of the diagram di and iii) 
replacing every loop generated in this process by n. The walled Brauer algebra BL,L{n) is 
the subalgebra of B2L{n) generated by the elements := £2^-1, Qs+i £2^, Pj,i+i := 
P2i-iP2zP2j-i, P^i+T P2iP2i+iP2j- The generators PjPi+iP^ are represented on the right 
in fig. [TJ For every diagram spanning B2l(ji) with vertices labeled as in fig. [TJ imagine 
moving all odd vertices to the left of a wall, all the even ones to the right, while keeping the 
connectivity unchanged. Then Bl^l^ti) is spanned by the set of B2L{n) diagrams, such that 
only strictly horizontal edge cross the wall. In this representation of BL L{n) we label the 
L up and L down vertices on the left of the wall by the set 1, . . . ,L from left to right and 
similarly those on the right by the set 1, . . . , Z. The P^i+i generators act on the left of the 
wall, the P^^TfY act on the right, while the generators Qn and Qu+i act across the wall. 

Next we give a brief description of a set of modules of BL,L{n) over which we actually 
numerically diagonalize the algebraic Hamiltonian H. These modules shall be related in the 
following to gl(n + A^jiV) traceless tensors of fixed co~ and contravariant shapes. For A 
and partitions of the non-negative integers |A|, = f < L, we denote by Ai,i(A,/i) the 
standard modules of Bj^ L^n). There are constructed in the following way. Let Sym(/) denote 
the symmetric group on / objects and S{X), S{^) denote its simple modules labeled by the 
corresponding partitions. Then, Al.l{^: m) has as basis the tensor products between the set 
of some diagrams on 2L points, with L of them on the either side of the wall, and some 
basis of S{X) and S{ii). The diagrams are such that every point is either free or belongs to 
a horizontal edge crossing the wall. The number of edges is fixed to L — f. We give a rough 
idea about how the action of the walled Brauer algebra can be constructed by diagrammatic 
multiplication in in fig. [5] i) in a first step, before the diagrammatic multiplication, assign 
labels to the free points of the diagrams on each side of the wall; ii) in a second step, after 
the diagrammatic multiplication, apply a surjective homomorphism from labeled diagrams 
to the tensor product of unlabeled diagrams with S{X) and S{p), that is to Al.l(A, /x). 
All labeled diagrams with more then L — f horizontal edges that can appear as a result of 
the diagrammatic multiplication belong to the kernel of this homomorphism. In particular, 
A/j(A,At) S{X) X S(^l). The detailed definitions can be found in |CVDM08| . 

We implement the action of BL,L{n) in Ai^i(A, /i) on a computer and investigate the 
spectrum of the algebraic Hamiltonian H. Before presenting and discussing the numerical 
results of sec. 13.31 one should explain how to extract from these data the spectrum of the 
original spin chain Hamiltonian H. 
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3.2 Traceless tensors 

Consider the gl(Af|iV) tensor {V'g)V*)^f. There is an obvious action of Syni(/) x Sym(/) 
on the V and V* factors. Let A, /i be partitions of /, which we symbohcaUy write as A, /i h /. 
One can apply a Young symmetrizer of shape A to the V factors and a Young symmetrizer 
of shape fi to the V* factors to get a tensor of shape t{X,fi) |Wey53| . We say that t{X,fj,) 
has rank (/,/), covariant shape A and contravariant shape fj, or, shortly, shape (A,/i). The 
number of Young symmetrizers of shape A is equal to the number of standard Young tableau 
of shape A, which is also equal to dimS'(A). Therefore, the number of tensors of shape 
(A, /i) is dim 5(A) x diniS'(/x). Every tensor t{X, /i) is a gl(M|A^) module, which appears as a 
direct summand in {V ® V*)^-^ . The symmetric group Sym(/) x Sym(/) acts in the space of 
tensors of the same shape (A, /x), transforming them into each other. The latter subspace is 
isomorphic to a direct sum of S{X) (8) S{^) modules of Sym(/) x Sym(/). These statements 
can be compactly written as follows 

(y^V*^^ dimi(A,Ai)5(A)®5(Ai) . 

Sym(/)xSyin(/) 

Consider now the subspace io(A,/i) C t(A,/i) of traceless tensors. Notice that this is 
a gl(M|iV) submodule, which will not necessarily be a direct summand of i(A,/i). More 
generally, one can consider the gl(M |iV) submodules f„_i(A, /i) C t(A. ^) composed of tensors 
whose all contractions of n covariant indices with n contravariant indices vanish. These 
provide a filtration of the tensor t(A, /i) 

to(A,M) C<i(A,//) c ••• Ct/(A,//) -i(A,/x) . (3.1) 

It is very important to observe that the subquotients <„(A,/i)/t„_i(A,/i) of this filtration 
are isomorphic to traceless tensors of lower rank (/ — n, / — n). For instance, taking a 
trace of ti(A,/i) one gets a traceless tensor of rank / — 1 because, by definition, all double 
contractions of ti (A, /i) must vanish. Now, the preimage of single traces of ti (A, fi) modulo the 
kernel to(A, fi) of the single trace homomorphisms is precisely to ti(A, /i)/to(A, /i). Therefore, 
ti{X, fi)/to{X, fi) is isomorphic to traceless tensors of rank {f — I, f — 1). 

We see that traceless tensors to{X,fi) have a fundamental role — all direct summands of 
the spin chain {V (E) V*)'^^ can be built out of them. As a consequence, the full spectrum of 
the spin chain Hamiltonian can be reconstructed from the spectra in the subspaces of traceless 
tensors of shape A, /i h /, / = 0, 1, . . . , L. In fact, not all of these shapes are possible. We 
shall determine the class of shapes for which the traceless tensors do not vanish later. 

The traceless tensors are not necessarily irreducible. One way to build submodules of a 
traceless tensor to{X, n) is by embedding into it quotients of traceless tensors of lower rank 
as follows. Let to{X' , fi') be a traceless tensor. A', ^' h / — fe. A' C A, /i' C and eA, denote 
some Young symmetrizers of shape X, fi. The tensor 

eAe^to(A', m') ® {{V ® y*)®^)^^*^!^) c t(A, ^i) 



10 



might have an intersection with a non-trivial submodule of to^X, fj,). The latter will generally 
be isomorphic to only a quotient of to{X',fi'), because the Young symmetrizers 6^,6^ are 
projectors. An illustrative example is the indecomposable gl(A''|iV) tensor = V 

V*. The traceless subspace to{l, 1), isomorphic to the adjoint representation, is spanned by 
elements of the form G^^ei^e^ subject to the constraint strG = G^(— l)'*' = 0, where {e^jf^]^ 
is a basis of V and {e*}f^]^ is the dual basis. The quotient l)/to(l, 1) is one dimensional. 
A coset representative for this quotient is, for instance, (— l)l'lei ® e*. The traceless tensor 
to(l, 1) is also indecomposable, but reducible. It has a unique proper non-trivial submodule 
spanned by the gl(A''|7V) traceless invariant Ci ® e'. 

Represent now the covariant part of every tensor t(){X, ii) of fixed shape (A,/i) and ranks 
(/, /) by / dots on the left of an imaginary wall and the contravariant part by / dots on the 
right of that wall. Then the walled Brauer algebra generators Pu+i will act on the left of 
the wall as in 5(A) and the P^^^ipj- generators will act on the right as in by transforming 
dim S{X) X dim S{iJ,) different traceless tensors of shape (A, fi) into each other. The generators 
Qii and Qu+i will act across the wall by contracting a covariant index with a contravariant 
one. In view of the tracelessness condition this action is trivial. Thus, the space of all traceless 
tensors io(A, /i) of the same shape (A, fi) is isomorphic to a direct sum of dim io(A, /i) modules 
A/j(A,/i) ~ S'(A) X S{fi) of the walled Brauer algebra. A very important observation is the 
triviality of the centralizer of a traceless tensor 

Endgi(MiAr) to{X, ^) ~ C . (3.2) 

This is a generalization of the Schur lemma for gl(n) traceless tensors, which are irreducible. 
The statement p.2p follows immediately from the action of the walled Brauer algebra in 
the space of traceless tensors of shape (A, /i) that we have just described. It means that 
traceless tensors are a very special type of indecomposables, namely any gl(M|A^) Casimir 
is diagonalizable and proportional to the identity in a traceless tensor. It should be noticed 
that this is typical of highest weight or Kac modules |Kac77b[|Kac77aj . 

Assumption 1 Traceless tensors are highest weight modules. 

This means that there is a Borel subalgebra fa of gl(Af |iV) ~ fa © and a fa-highest weight 
vector V e to(A,/i) such that the full tensor to{X, fi) can be generated from v by repeated 
action of n^. We shall see later how to choose fa for given to (A, /i). 

Consider now the vector space 6l,l{X, fj) of all possible embeddings of traceless tensors of 
shape (A, fi) and ranks (/, /) into the spin chain {V ® V*)'^^ . It consists of tensor products 
of tensors to (A, /i) with gl(M|iV)-invariants of {V ®V*)®^^~f\ Notice that there is a unique 
gl(M|A^) invariant in V ®V* , which can be written as Ci® . Representing this invariant 
by an edge connecting two vertices across the wall and the traceless tensors to(A,/i) as we 
did before, we can visualize 5l,l{X, fJ-) as a diagram with L vertices on each side of the wall 
and L — f edges connecting pairs of vertices across the wall. Thus, we reconstruct the same 
diagrammatic representation of (5l^l(A, fj,) as for Al l(X, fi). This proves that all the relations 
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between the generators of B^xin) satisfied in Ai^i(A, /i) will be satisfied in 6L,L{^,fJ-) as 
well. The converse is generally not true, meaning that SL,L{^,fJ') is generally only a quotient 
of {dimto{X, fi) direct sums of) Ai^i(A,/i). We stress that neither 5l,l{X,Ij) nor Ai^i(A,^) 
are necessarily simple BL^iji) modules and, therefore, the quotient can be non-trivial. 

The ABA in some grading E provides gl(A/| A^) Bethe eigenvectors of highest weight with 
respect to the Borel subalgebra by, determined by E and the ordering (12.71) . Therefore, in 
order to match the numerical spectrum of H in AL.L(A,/i) with the exact spectrum of H by 
the ABA we need to know the highest weight of a traceless tensor to{X, ij,) at least in one 
grading E. We explain below how to evaluate it. 

Consider the Young diagrams corresponding to the shape {X, fJ-) of a full tensor t(A, /i). 
The basis vectors in the tensor subspace of co(ntra)variant shape A (/i) can be represented by 
co(ntra) variant Young supertableaux of shape A (fi), that is Young diagrams of shape A (p) 
with a fundamental weight (— Cj) inscribed in every box. The pattern of weights within the 
Young diagrams must satisfy the supersymmetrization rules, that is i) in the same row the 
bosonic (fermionic) weights are weakly (strongly) ordered w.r.t. each other, ii) in the same 
column the bosonic (fermionic) weights are strongly (weakly) ordered w.r.t. each other and 
iii) bosonic weights are weakly (strongly) ordered w.r.t. fermionic weight in the same row 
(column). The weight of a supertableau is the sum of all weights it carries in its boxes. 

Both A and fi must fit into a hook whose horizontal (vertical) arm is of width M (N) , 
otherwise i(A, fj.) vanishes identically because it is not possible to fill in the Young diagrams 
and get Young supertableaux compatible with the supersymmetrization rules. The highest 
weight of a supertableau depends on the grading. The choice of grading is a splitting of the 
set of basis vectors into two sets B [F] — {e^ | \i\ = 0(1)}^ which are ordered w.r.t. the 
total ordering (j2.7p . Equivalently, it can be represented by paths connecting the two corners 
of the (M, Af )-hooks as represented in fig. |31 Fix these paths and consider Young diagrams 
A, /i with rows Ai, /ii and columns A^, /i^. Then, the highest weight of t(A,/i) can be written 
in the following form 

M N 

Ae(A, h) = ^heb(,) - f,e^,)] + ^[Qe/(,) - Qej(,)] (3.3) 

i=l i=l 

where h{i) and f{i) are the elements at position i in the ordered sets B and F, b{i) and f{i) 
are the elements at position i in the ordered sets ~B and — F, — max(0, Xi — X]j=i(l ~ 
(-l)l^l)/2), c, = max(0, A^ - Y.U{1 - (-l)l'l)/2, r" = max(0, - Y^^il - (-l)l^l)/2) 
and Ci = max(0, ^'^ — + (^l)'*')/2) are number of boxes in a row or column of A 

or /i overpassing the grading paths as represented in fig. [31 

The highest weight component of a tensor i(A,/i) w.r.t. the grading E will belong to 
the traceless subspace io(^7 A*) if the highest weight Young supertableau of shape (A, /i) does 
not contain some fundamental weight ±6^ in both A and /i. Otherwise, the corresponding 
highest weight component of i(A, /i) w.r.t. the grading E will either i) not belong to to(A,/x) 
or ii) generate a submodule of tQ{X, fi) isomorphic to the embedding of a quotient of a lower 
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Figure 3: Covariant Young tableau A and contravariant Young tableau /i of (5, 4)-liook shape 
for gl(5|4). Bosonic fundamental weights belong to B — {e2, €3, 65, eg, eg}, while fermionic 
weights to F = {ei, £4, ej, eg}. Correspondingly —B = {—eg, —eg, — e5, —£3, —£2} and —F = 
{-£8, -£7, -e4, -ei}. The values of r4, rs, C3, C4, r2, • ■ • jT^s, C3, C4 are zero. 



rank tensor to{X',n')- If the latter case holds, then the possible Young diagrams (A',^') 
are obtained from the Young diagrams (A, fi) by removing pairs of boxes from the highest 
weight Young supertableau of shape (A, /i): a box of A carrying some weights £i together with 
a box of /i carrying the opposite weight —ei. Moreover, in the case ii) the highest weight 
of the top0 of io(A,/i) will be smaller then the highest weight of the submodule to(A',/i')- 
Therefore, io(A,/i) cannot be a Kac module w.r.t. bs- 

We call g\{M\N)~-admissible the shapes (A,/i) for which the gl(M|iV) traceless tensors 
to{X, fi) neither vanish nor are isomorphic to lower rank traceless tensors. We are now ready 
to answer the very important question: what are the admissible shapes of traceless tensors? 
According to the previous discussion on the highest weight component of a tensor i(A,/i), 
a shape (A, /i) is admissible if there exists a grading S such that the highest weight Young 
supertableau of shape (A, /i) does not contain any fundamental weight ±£j both in A and in 
11. Consequently, for a traceless tensor to{X,fi) of admissible shape there is a grading T, and 
a corresponding highest weight (j3.3D such that no cancellation between the and or Ci 
and Ci terms can occur. If ?Js(A), denote the number of nonzero "reduced rows" r^, fj 

and /is (A), h^{^) denote the number of non-zero "reduced" columns Ci, then one must 
have 

3Y. such that we(A) + < M, hY,(\) + hT.{^J) < N (3.4) 

for a gl(MlA^)-admissible shape (A, /i), as represented in fig.[3]These admissible shapes nicely 
reduce to hook shapes for purely covariant or contravariant gl(A/| A^) tensors and to staircases 

^The top of a module is the quotient by the intersection of all maximal ideals. The top of a Kac module 
is the irreducible representation of highest weight. 
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for gl(n) traceless tensors [BCH+qI) . 

We say that a shape (A, /i) of a traceless tensor to{X, fj.) is E-admissible if the inequahties 
in eq. p.4p are satisfied. Let K^{A) be the Kac module of highest weight A w.r.t. bs. We 
can make now assumption 13.21 more precise. 

Assumption 1' The following isomorphism holds for S-admissible shapes (A, /i) 

to{X, fi) - K^{A^{X, fi)) . 

This assumption in combination with the general theory of Kac modules jKac77b[|Kac77aj 
is very useful for counting highest weight vectors. Namely, if (A,/i) is S-admissible then the 
number of highest weight vectors in to{X, fi) w.r.t. bs is equal to the number of irreducible 
subquotients. Noticing that a highest weight vector cannot generate more then a highest 
weight module, we prove the following claim in app. |^ 

Claim 1. Highest weight vectors of {V ® V*)^^ w.r.t. any Borel subalgebra belong to sub- 
modules isomorphic to traceless tensors to{X, fi), A, /i h / = 0, 1, . . . , i of admissible shape. 

To sum up, we have explained the connexion between traceless tensors to (A, /i) of admis- 
sible shapes (A, /x) and standard modules Al^l{X, /i). Secondly, if (A, fj.) is S-admissible, then 
eq. (13.31) allows to compute the highest weight of to{X, fj.) w.r.t. the Borel subalgebra bs. 
Evaluating the highest weight of tQ{X, ii) w.r.t. arbitrary Borel subalgebras is more delicate, 
mainly because of indecomposability issues. Finally, claim [T] indicates a natural relation- 
ship between traceless tensors and Bethe vectors, which have the highest weight property, 
constructed in the framework of ABA. 

3.3 Spectrum 

We present the spectrum of H in various Al,l(A, /i) standard modules of BL^Lin) in tab. 13.31 
At a first glance, it appears that the vacuum always lies in Ai^i(0, 0). To check more thor- 
oughly this vacuum hypothesis we need an additional assumption on the spectrum. Consider 
the spectral sets of H defined as 

spec/= y specAL,L(A,/z) . 

Notice from tab. 13.31 that the lowest eigenvalue in spec / always lies in Al.l(1''', I'''), where 
1^ denotes the Young diagram with a single column of length /. We have checked this 
observation extensively (more details). 

Assumption 2 The lowest eigenvalues of H in spec / always lies in , 1^). 

Comparing only the lowest eigenvalues in Al.l(1''^, 1'^) allows us to gain several spin chain 
length units and check the vacuum hypothesis further, see tab. 13.31 

To extract the spectrum of the gl(n -I- N\N) spin chain Hamiltonian (I2.20p from the 
spectrum of the algebraic Hamiltonian H we do the following. For every pair of Young 



14 



Table 1: Lowest eigenvalue of H in Al(A, with opposite sign. 



L 


n 


(P,2) 


(2,2) 


(1^21) 




(21,21) 


(21,3) 


(3,3) 




1 


32.130 


29.797 


29.461 


27.453 


27.156 


24.101 


21.446 






±0.591i 
















2 


23.237 


21.904 


22.626 


20.498 


20.506 


17.999 


15.548 


r 




±0.339i 














5 


3 


20.763 


20.130 


20.554 


18.530 


18.762 


16.572 


14.472 




4 


19.968 


19.487 


19.588 


17.682 


18.139 


16.180 


14.153 




1 


41.558 


39.639 


38.640 


36.835 


36.559 


33.762 


31.264 










±G.355i 




±0.513i 


±0.127i 






2 


29.537 


28.282 


28.802 


26.976 


26.984 


24.784 


22.607 










±G.267i 




±0.426i 


±0.148i 




6 


3 


26.123 


25.600 


25.884 


24.210 


24.365 




20.713 










±0.194i 




±0.258i 


22.550 






4 


24.917 


24.557 


24.566 


23.031 


23.416 


21.800 


20.021 










±0.129i 











diagrams (A, /x) with / = 0,1,..., L boxes and of admissible shape we pick a grading S 
such that the inequalities in eq. (13.41) are satisfied. Then we try to reproduce the spectrum 
of H in Az,_i(A,^) by means of eq. (|2.18p from numerical solutions of BAE (j2.10l) in the 
the form determined by the grading E and for root numbers corresponding to the highest 
weight As(A,^) in eqs. (|2.121 [H^]) . rl Not all eigenvalues can be reproduced in this way. This 
happens because, as we have explained in sec. 13.21 the vector space 5l,l{^tI^) of all possible 
embeddings of gl(n + N\N) traceless tensors to(A, /i) into {V ® 1/*)®^ can be identified with 
only a quotient of the standard module Al.l(A, ^) of BL^L{n). When an eigenvalue of H can 
be reproduced this way, we assume that a corresponding non-vanishing Bethe vector exist. 
Otherwise, we assume that there is no eigenstate of H corresponding to that eigenvalue. So, 
we rely on the completeness of the ABA, at least as far as the spectrum is concerned. 

Denote the spectrum of the gl(M|A^) spin chain Hamiltonian ()2.20p by s^igc H m\n ■ The 
central idea of this section was the existence of an abstract algebra BL j^{n) such that the cen- 
tralizers of the series = 0, 1, 2, . . . of gl(n + A^|A^) spin chains C{L) provide different, A^ de- 
pendent, representations of ^^^^(n). This suggests that the intersection n7vez+ spec Hn+N\N 
might be non-trivial. In fact, with the cohomological techniques developed in |CCMS10] one 
can prove the following relationship between the spectral sets spec iJ^_|_„|7v with n fixed 

speciJ„|o C speci?„+i|i c speciJ„+2|2 C • • • C H . (3.5) 

This "embedding of spectra" is a very interesting and general feature of supergroup spin 
chains and one might wonder how does it carry on to the field theory description of the 

^One could work with a single form of BAE, say, that corresponding to the distinguished gradation Sq. 
However, the weight of Bethe vectors reproducing eigenvalues of H in Ai^^l{\, fi) will no longer be given by 
eq. I I3.3I I if the shapes of (A,/i) do not satisfy the inequalities 113.411 w.r.t. Eq. 
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Table 2: Lowest eigenvalue of H 


in Ai(l''', l'^') with opposite sign. 




L 


n 


fc = 


k = 1 


fc = 2 


/fc = 3 


k = A 


k^b 


k^Q 


fc = 7 


fc = 8 




1 


40.000 


38.846 


36.564 


32.710 


26.385 


20.000 










2 


28.062 


26.369 


26.156 


25.082 


22.606 


20.000 








5 


3 


24.625 


22.950 


22.944 


22.640 


21.407 


20.000 










4 


23.123 


21.631 


21.456 


21.456 


20.828 


20.000 










1 


48.000 


46.723 


45.626 


41.538 


36.782 


30.397 


24.000 






6 


2 


33.550 


32.126 


32.126 


30.902 


29.253 


26.647 


24.000 






3 


29.388 


27.989 


28.024 


27.617 


26.885 


25.476 


24.000 








4 


27.574 


26.339 


26.146 


26.086 


25.754 


24.921 


24.000 








1 


56.000 


55.134 


53.631 


50.894 


46.015 


40.796 


34.399 


28.000 




7 


2 


39.054 


37.826 


37.699 


36.992 


35.260 


33.307 


30.660 


28.000 




3 


34.172 


32.971 


32.939 


32.748 


31.954 


30.981 


29.504 


28.000 






4 


32.046 


30.991 


30.800 


30.800 


30.420 


29.885 


28.962 


28.000 






1 


64.000 


63.035 


63.296 


59.127 


55.589 


50.296 


44.799 


38.400 


32.000 




2 


44.569 


43.488 


43.488 


42.533 


41.457 


39.453 


37.325 


34.664 


32.000 


8 


3 


38.970 


37.917 


39.480 


37.652 


37.146 


36.134 


35.019 


33.515 


32.000 




4 


36.530 


35.611 


35.428 


35.391 


35.169 


34.603 


33.944 


32.982 


32.000 



continuum limit. In this respect, two scenarios are possible. The first possibility is that 
spec i7„_^^/|jv/ becomes a very excited subset within spec_ff„^jv"|JV": where N' < N" , and 
decouples from it in the thermodynamic limit i — s> cxd. This means that the vacuum energy 
per site of Hn+N'\N' is higher then the vacuum energy per site of Hn+N"\N" in the thermo- 
dynamic limit. The second, more interesting, possibility is that the vacuum energies per site 
for both Hamiltonians coincide in the thermodynamic limit. The first possibility occurs, for 
instance, in the V®^ chains with Hamiltonian ±J2Pi,i+i- 

In next section we describe the mechanism which is behind the embedding of spectra p.5p 
at the level of BAE. The answer to the question of how excited spec i?„_(_jv'|7V' is within 
spec i7„_(_^//|7V" for N' < N" will have to wait until sec.|H 

4 Restriction and lift of BAE 

At the end of the previous section we have explained the phenomenon p.5p of embedding 
of spectra. For it to work, it was important that there is a series {gl(n + N\N)} of 
spin chains endowed with Hamiltonians Hn+N[N which are all different representations of 
the same algebraic Hamiltonian H e -Bl,l(?t-)-u The integrability did not matter. In this 
section we would like to understand how this curious phenomenon arises at the level of BAE 
describing the spectra of very general gl(M|A^) integrable Hamiltonians and what structure 
is responsible for it. 

Fix a grading = {a, = (-l)l*l}fl|^ of the ordered basis vectors ^TT^ of the gl(M|A^) 

^Strictly speaking, it does not make sense to talk about different representations of a single matrix H. 
What we mean here is that the i?„-)-jv|JV ^'"^ the image of H in different representations of Bi^j^in). 
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fundamental representation V and let A^'''^ — {aj = ej — be the corresponding 

simple root system, where r — M+N ~1 is the rank. We have explained how to construct the 
ABA in the grading E for the spin chain [V 'iSiV*)'^^ in secjH The Bethe vectors are highest 
weight vectors of weight (|2.12l) w.r.t. the Borel subalgebra bs determined by the grading E 
and the ordering (j2.7p . that is by Aq^'^. We shall restrict to Bethe vectors whose weights 
are given by highest weight Young supertableaux of E-admissible shapes (A, /i) according to 
eq. p.3p . Due to fundamental role of traceless tensors explained in sec. 13.21 and claim (TJ it 
is clear that with the imposed restriction, one must consider the ABA in all the gradings 
in order to recover the full spectrum of the spin chain Hamiltonian H . For a fixed grading 
and a corresponding simple root system, fig. [3] and eq. (j3.3l) implies that we are restricting 
to Bethe vectors of weight w = X^t^i'^ w^^i, such that for every simple bosonic root 

aj = ej — ej+i one has > w^^^ , while for every simple fermionic root one has the following 
implications 

< — > w^+^ < , =0 — > < (4.1) 

w^+^ > — > , w^+'^ =0 — >w^ > . 

The main reason for introducing this restrictions and working with BAE in multiple gradings 
is the bounds on the number of Bethe roots resulting from eqs. (14. ip . 

4.1 Restriction 

We wish to consider the BAE in the form (|2.9p corresponding to a simple root ak such that 
the following assumptions hold 

Al ak is odd, that is akCk+i = —1 

A2 ak has no source terms or, equivalently, k ^ l,r 

With these assumptions, we have v^''^ BAE for of the form 

ak-iA.k-i{uf^^^) ^ Qk-2{uf^^'')Qk-i{uf^^^ + (Jk)Qk{uf'^^ -o-fc) 



akAkiuf-'^) Qk-2{uf-'^ + ak-i)Qk-iiuf~'^ - CTk-i)Qk{uf-'^) 



v'^^^ equations for ak 



and i/C^'+i) equations for 



Qk-iiuf^)Qk+i{uf^ +ak) 

Qk-i{uf + (Tk)Qk+l{uf) 



Otk+l 



akAk+duf+''>) _ Qk{uf+'^)Qk+i{uf+'^ + ak+2)Qk+2{uf+'^ - cTk+2) 



Uk+2Ak+2\u) ') Qk[Uj ' - (Tk)Qk+l(u) +(Tk)Qk+2[u) ') 



(4.2) 



(4.3) 



(4.4) 



Notice that if j/C^-i) = jyC^'+i) then one can reduce the BAE (j2.9p for the gl(M|A'^) spin chain 
{V®V*)®^ to the BAE for the gl(M - 1|7V - 1) spin chain of the same type {V ®V*)®^ by 
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Rl identifying the Bethe roots corresponding to ak^i and ak+i 

wf-''}^^r - {-r^^i^r' (4.5) 

R2 multiplying the BAE for ak~i and ak+i corresponding to, say m^'^ = '"j*^^^'' 

Crfc+2Afc+2("j ) Qk-2[Uj + (Jk-l)Qk-l{u) - (Jk-l)Qk+2[Uj ') 

(4.6) 

Indeed, eq. (|4.3p is trivially satisfied because according to Rl one has Qk-i{u) — Qk+i(u). 
Furthermore, multiplying BAE according to R2, the Bethe roots {itj'''*}jli drop off yielding 
a BAE of the form (|4.2p with k replaced by A: + 2. We call a restriction of BAE the procedure 
Rl R2. 

The restriction procedure can be given an algebraic meaning and, therefore, partially 
explained as follows. As described in detail in jCCMSlO] . choose 

Q = Ek+ik (4.7) 

to be the odd gl(M|A^) element that squares to zero and defines the gl(M — 1|A^ — 1) chain 
(VAf-i|Ar-i ® V7f_i|^_i)®^ as the Q-cohomology of the g\{M\N) chain {Vm\n «> V^j^j^)'^^. 
A necessary condition for a highest weight vector uj £ {Vm\n ®Vm\n)^^ *° yield a non-trivial 
Q-cohomology is for it to be in the kernel of Q and, therefore, one must have 

[Ekk+i,Ek+ik](^ = (wt(w),afe) = . (4.8) 

Applying this constraint to a Bethe vector of weight (j2.12D one recovers the condition z/C^"^) = 
necessary for Rl to hold. The form of the reduced gl(M — 1|A^ — 1) BAE can also 
be easily understood. The subquotient (ker^^^^ Ce^) /Cafc of the gl(Af|A^) weight 

space can be straightforwardly identified with the gl(Af — 1 1 A^ — 1) weight space J2i^t^+i Ccj- 
Therefore, the subquotient ( ker^^^.^^) Aq^'^) /Cak of the gl(Af | A^) simple root system Aq^'^ 
induced by cohomological reduction can be identified with a gl(Af — 1|A^ — 1) simple root 
system Aq^ ^'^ ^ = {ai, . . . , afc_2, efc-i — efc+2, afc+2, • ■ • , Q^r}- The reduced BAE have the 
form (|2.10p corresponding to precisely the simple root system Aq^ ^'^ ^ . 

4.2 Lift 

To resume, assuming A1-A2 holds for the BAE of the spin chain {Vm\n '55 ^m\n)'^^ written 
w.r.t. a simple root system Aq^'^, we showed that one can restrict them to the system of 
BAE for the spin chain {Vm~i\n-i ® ^m-i\n-i)'^^ written w.r.t. the simple root system 
Aq^ i|JV 1 _ 1^^^ ^ afc_2, Efc-i — efc+2, ak+2, ■ ■ ■ , otr} induced by cohomological reduction. 
This restriction is represented at the level of Dynkin diagrams in fig. |4l 

(k) 

Notice that we have not imposed any condition on the roots Uj . Therefore, it is le- 
gitimate to ask if the BAE satisfying A1-A2 actually admit solutions of type Rl. So, we 
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ak-1 at Ofc+i 

Figure 4: Restriction of BAE for supergroups. The dashed hues indicate possible roots on 
the left or on the right. For a root aj = ej — tj^i with j = k,k ±1, we have indicated with 
± signs the gradings of fundamental weight determining it. 

are given a simple root system Aq^'^ and a solution / = 1, . . . , fc — 1, fc + 2, . . . , r 

of the {Vm-i\n-i ® ^M-i\N~i)^^ BAE w.r.t the simple root system Ag^^"^'^^"^ induced 
by cohomological reduction from Ag^'^ with Q as in eq. (|4.7p . We must show that there 
is a solution {wf jj^i, ? = 1, . . . , r of the {Vm\n ® V^mi^v)®"^ BAE w.r.t. A^'^ which re- 
stricts to the given one. First of all, eq. (j4.8p implies i/C^+i) = jyC^-i) — p. After defining 
M^'^^^-' = uj'' for j = 1, . . . , P the task is reduced to finding a solution to either 

eq. (|4.2p or eq. (|4.4p . where we have set p = It p > P, such a solution obviously exists. 
More then that, for p > P there is a continuum of such solutions! However, recall that for a 
fixed root system and corresponding BAE we have restricted to Bethe vectors such that their 
weights satisfy the constraints (14. ip . If the solution we are looking for exists then the weight 
of the corresponding Bethe vector lu can be written as wt(a;) = • • • + {P — p)(ek — Cfc+i) + . . . . 
The constraints (j4.ip imply p < P. So, for p — P the solution always exists. In our numer- 
ical investigations we have observed that solutions might exist even for p < P. We call the 
solutions of g\{M\N) BAE constructed in this way from solutions of g\{M — 1|A^ — 1) BAE 
lifted solutions. We discard lifted solutions with p < P for the following reason. 

Claim 2. // a solution of g\{M - l\N - 1) BAE admits lifts to solutions of g\{M\N) BAE 
corresponding to Bethe vectors of different weights, then all of them have zero Q-cohomology. 

We prove the claim in appendix |B] 

Next, we show that the lift is unique. So, we have seen that lifted Bethe vectors correspond 
to lifted solutions with v^''^ — zy('^^i) — P. Let us now prove that the P equations, say, (|4.2p 
regarded as a constraint on the unknowns admit a unique solution. We start by 

multiplying both sides with Qk{u^j'~^^) and then expand in {u^'^^jj^^. What we get is a 
system of P linear equations for the P unknowns 

E(fc) (k) 

il<-<il 

Thus, for a given solution {uf^}fl\, j = 1, . . . , fc - 1, fc 2, . . . , r of the gl(M - l|iV - 1) 



gl(M- l|Ar- 1) 



p 




-0-- 
ttfc^i ~ ttfe+i 
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BAE, this system of linear equations admits a single solution for the {sjl^Li- Finally, notice 
that the latter determine the set {Wj'^''}jLi uniquely, because every Bethe root in this set is 
a solution of the polynomial equation of degree P 

p 

Y[{u - uf) = + + ■ • • - {-ifsp-iu + {-ifsp = . 

There are important cases when one can compute the roots explicitly in terms 

of the other roots. This happens if 

C 1.1) ak±i is odd and sourceless 

C 1.2) lyC'^^') = P, 

where we meant that either we choose the plus sign or the minus sign overall. In the case 
C 1 eq. (|4.2p or eq. (|4.4I) can be rewritten in terms of variables (I2.11[) as 



^(fc±l) _ ^(fc) ) = Jl p- fjk±l} _ Jk±2) ^ 
i=l i=l 

We immediately read off the obvious solution 

{4'^}f.i-{xf^^'}f.i, (4.9) 
which we already know is unique. Another important case is 
C 2.1) fc = 2, ai is odd and P = L 
C 2.2) fc = r - 1, is odd and P = L. 
In the case C 2.1 the one has 

i—l a—1 

while in the case C 2.2 

L L 



where we have used eqs. (j2.10p with non-zero arbitrary inhomogeneities. Again, we read off 
the unique solutions 

W'hf^i = {ya}Li , {4~'^}f=i = {y-a}a=i ■ (4.10) 

It is important to notice that inhomogeneities are essential to lift the degeneracy of solutions 
in the case C 2. 
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4.3 Generalizations 

We have argued that the restriction and hft of BAE have an algebraic origin. From this 
viewpoint, assumption A2 seems unnecessary. Indeed, all of the above constructions can be 
appropriately modified to accommodate the boundary case corresponding to Al and 

A2' fc = l,r. 

The necessary condition (|4.8p for a non-vanishing Q-cohomology of a gl(M|iV) Bethe vec- 
tor uj implies that reducible solution of g\{M\N) BAE must have j/^^^ = L if fc = 1 and 
j^(r-i) — 2v if fc = r. The reduced BAE must have a form (|2.10l) corresponding to simple root 
systems Aq^ ^'^ ^ — {013, . . . , a^} and A^"^ ^'^ ^ — {ai, . . . , ar-2} induced by evaluating 
the subquotient (ker^^^^ _^ Aq^'^)/Cq!A; for fc = 1 and k = r respectively. To satisfy these 
requirements one must replace Rl— R2 with 

R' = {2/a}Li for k = 1 and {x^'-'^}f^, = {yalLi for k - r. 

Invoking again claim [5] and the arguments that follow it, we conclude that the lift of a 
gl(M — 1|A^ — 1) Bethe vector to a gl(M|A^) Bethe vector must be unique and satisfies 
ly^^^ = L for k = 1 and j/^''^ = L for k = r. Finally, the corresponding lifted solution can be 
evaluated explicitly if a2 is odd for fc = 1 and a^-i is odd for k = r 

Subsectors in integrable systems are not a new phenomenon. For example the homoge- 
neous gl(27i + 1) spin chain {V (E) V*)'" of sec. [2] contains a subsector corresponding to the 
2L-th tensor power of the fundamental representation of osp(l|2n), see |SWK02j . We shall 
provide more examples in sec. [5] However, all these examples lack by far the generality of the 
restriction and lifting procedures we have described. This is because in our case there is an 
algebraic mechanism behind which relates the representation theories of {gl{n + N\N)}jffzi+ 
Lie superalgebras. The same cohomological reduction mechanism exists for any Lie super- 
algebra except osp(l|2ri). Therefore there is no doubt that the same restriction and lift 
phenomena occur also for general osp(i?|2S') integrable spin chains. Moreover, the embed- 
ding of spectra p.Sp was already observed on a (non-integrable) osp(2S' -I- 2125*) spin chain 
in |CS09a| . 

We come back now to the spin chain of sec. [5] and answer the question about how excited 
w.r.t. each other are the subsectors (13.51) . 

5 Vacuum and low lying excitations 

In this section we present numerical evidence showing that all integrable gl(n -I- A^jA'^) spin 
chains (|2.20p with n fixed have the same vacuum energy. In view of the embedding of 
spectra p.Sp . we introduce the notion of degree of an excitation, which is the smallest value 
of N for which it appears in spec if„+jv|Ar- We then proceed to classify the excitations of 
degree and 1, and present the form of the numerical solutions of BAE reproducing them. 
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5.1 Bosonic lift 



From (j3.5p , the vacuum energy of the spin chain is greater or equal to the vacuum energy of 
the algebraic Hamiltonian. Supposing they coincide, from tables p. 31 13.31) and the general 
discussion of sec. [3] it follows that the vacuum state in the spin chain is a gl(n + N\N) 
invariant tensor. According to eq. (|2.12p . the number of roots in the vacuum state should 
then be i/*-'^^ = L. Numerically, we confirm that there is indeed a solution to the gl{n + N\N) 
BAE (I2.10p with i^''^^ = L which reproduces the vacuum energy of the algebraic Hamiltonian 
H G BLj^{n). According to the terminology of sec. HI this is the gl(n + N\N) lift of the 
vacuum solution for the purely "bosonic" gl(n) chain. The latter is characterized by n — 1 
seas of real roots without holes. 

The low rank cases gl(2|l) and gl(3|l) have to be treated independently. For the gl(2|l) 
case it was realized in [EFS05I that the vacuum energy of the Hamiltonian (|2.20p is exactly 
—AL and the corresponding vacuum solution in the grading S = {+,—,+} is highly degen- 
erated x^^^ = xj^'' — 0. This solution obviously does not make sense in the framework of the 
ABA. The problem persists in all other gradings as well. A regularized vacuum Bethe vector 
can be defined by introducing arbitrary inhomogeneities for every V and V* site in the chain. 
Then, according to eqs. (I4.10|) . the vacuum solution is entirely fixed by these inhomogeneities 
and the vacuum Bethe vector can be constructed. The vacuum state of the homogeneous 
chain is then defined by a limiting procedure. The gl(3|l) spin chain has the same problem, 
again in all the gradings, and the same cure. 

We now realize the embedding speciJ„|o C spec _ff„+jv|Af explicitly. Following [CCMSIO] . 
choose an odd element Q G g\{n + N\N) that squares to zero, has rank N and defines the 
gl(n) spin chain V*)®^ as the Q-cohoniology of the gl(n-|-iV|7V) spin chain {V^V*)'^^. 
For this choice of Q, one can perform the lift of all gl(n) Bethe vectors to gl(n + iV| A^) Bethe 
vectors as described in sec. E) As we shall see in a moment, the lift depends on the choice 
of Q. We call this part of the spectrum the bosonic lift. The latter is straightforward to 
understand in terms of gl{n) spin chain excitations |AL80a[rAL80bl[AL79[|DN98j . 

How to characterize the gl(n + N\N) symmetry of the excitations that belong to the 
bosonic lift? We say that the component Wi of a gl(n + N\N) weight w = X^iH^^^*^* 
bosonic if |i| = and fermionic if \i\ = 1. According to sec. 21 if a gl(n + N\N) Bethe 
vector u! is the lift of a g\{n) Bethe vector, that is uj has degree 0, then it should be possible 
(in some grading) to write its weight in the form wt{uj) = '^fLt^ Witi with no fermionic 
components and at most n bosonic components. Equivalently, the weight of a Bethe vector uj 
of degree must be representable (in some grading) by highest weight Young supertableaux 
of gl(n)-admissible shape (A, /i) 



These shapes are E-admissible and most obviously gl(n) reducible w.r.t. the gradings 



(5.1) 




2N 



J ; 



(5.2) 
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Figure 5: The grading [+'"(H — )^+"^™], 1 < m < n — 1 is such that there is a chain of 2N 
odd roots. On the left: lift of the gl(»T.) vacuum to the gl(?T. + N\N) vacuum. On the right: 
Bethc roots for Bethe vectors in the bosonic lift w.r.t. the two choices (15.31) of Q. 

where \'i < m < n — fi^. There are many choices for the nilpotent element to define the 
restriction and lift in the grading (j5.2p . For instance, the two possibilities 

Ql E CE„i+2,m+l © CE,n+4,m+3 ® • ' " ® 'CEjn+2N ,m+2N ~1 (5.3) 
Q2 G CE,n+3,m+2 © Cii^„i+5,m+4 © ' • • © CEm+2N+l,m+2N 

are represented in fig. [5l Notice that for a fixed Q the lift is unique, although it is different 
for different Q's. Multiple gradings of the type (|5.2p must be considered in order to recover 
the full bosonic lift. 

To conclude, we write down explicitly the vacuum solution in the grading [+™(H — 
with 1 < m < n — 1 

4") = 4'"+^' = • • . = 4™+^^) (5.4) 

_ ^('"+3) _ _ (m+2W+l) 

where the roots {x^'''}^^!, I — \, . . . ,m,m + 2N + 1, . . . , n + 2N — 1 correspond to the gl(n) 
vacuum solution. We represent the lift of the bosonic vacuum in fig. [5] We have already 
mentioned, see eq. (|4.10|) . how to treat the special cases gl(2|l) and gl(3|l) by introducing 
arbitrary inhomogeneities. 

5.2 Degree one excitations 
5.2.1 Classification 

Fix an odd element Q G gl{n + N\N) that squares to zero, has rank — 1 and defines 
the gl(n + 1|1) spin chain {V © T/*)®^ as the Q-cohomology of the gl(n + N\N) spin chain 
(V © V'*)®^. Excitations of de gree 1 correspond to gl(n + N\N) Bethe vectors which: i) 
are non-reducible to gl(n) Bethe vectors and u) are lifted gl(n + 1|1) Bethe vectors. We 
shall call the corresponding Bethe vectors also of degree 1. According to sec. SI the weight 
of a Bethe vector oj of degree 1 admits a representation (in some grading) of the form 
wt(a;) = '^iLt^ WiCi with at most rt + 1 bosonic components and at most one fermionic 
component. Notice that due to the claim [2] of sec. |4l there are also excitations of degree 1 
with the same weight as excitations of degree 0, that is no fermionic components and at most 
n bosonic components. Taking into account sec. 13.21 as well, the weight of Bethe vectors of 
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Figure 6: Root numbers for gl(n + iV| A^) excitations of degree 1 and shape (A, /i) satisfying, 
top left, eq. w.r.t. some grading with < m < n + 1 - X'^ ot, 

bottom left, cq. dill) w.r.t. some grading [(+)'" (+)""'"] with X[ < m < n + I ~ fi'^. 
On the right: root numbers for the lowest lying degree 1 excitation. 

degree 1 must be representable (in some grading) by highest weight Young supertableaux of 
gl(ri + l|l)-admissible shapes (A, /i), that is at least one of the two conditions holds 

A'l + /i2 < " + 1 (5.5) 
A'a +n[<n+l . (5.6) 

Shapes (15. 5p are S-admissible and most obviously gl(ri + 1|1) reducible w.r.t. the gradings 
[(+)"-'"(+-)2^(+)'"], such that fj.[<m< n+l-A^, and the choice Q e CEm-n+3^m-n+2® 
■ ■ ■ (BCEm-ii+2N~i,rn-7i+2N~2- Shapes (15. 6|) are S-admissible and most obviously gl(n + l|l) 
reducible w.r.t. the gradings [(+)"'( — h)^^(+)""™], such that A'^ < m < n + 1 - ^2: and the 
choice Q G CEm+3.rn+2® ■ ' ' ® E„i+2N -i,m+2N -2- The Corresponding root configurations are 
represented in fig. [5] on the left. 

5.2.2 Low lying excitations 



Assumption 2 of sec. 13.31 and tab. 13.31 clearly indicates that the lowest lying excitation of 
the BL^Lin) algebraic Hamiltonian which is of degree higher then lies in Ai^i(l'^, l'^) 
with fc = [^] + 1. Equivalently, the lowest lying excitation of the gl(7i + A^|A^) spin chain 
Hamiltonian which is of degree higher then is a traceless tensor ^0(1*^, 1*^)- Eqs. (|5.5[ 15. 6p 
imply this excitation is of degree 1. To reproduce it from gl(rt + 1|1) BAE wc have chosen 
the following gradings 

„ f , n = 2m-l 

11 o (^-^^ 

I n ~ Zm 



The corresponding root numbers are represented in fig. [5] on the right. We have performed 
extensive numerical calculations on the form of the solution of Bethe ansatz equations re- 
producing the lowest lying excitations of degree 1. Compared to usual Lie algebra spin 
chains, this excitation looks very strange, to say the least. Strikingly, for n odd the roots 
{x'-^^}f-{^, {x5™+^^}j^rr are always complex, while for n even the roots {x'f'^^^}f-^-'^ 



^We will show a plot of the root distributions in the complex place for low lying excitations of degree 1 in 
a future version of this eprint. 
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are always complex. We recall that all roots are real for the vacuum solution (I5.4p . Thus, 
this excitation is not constructed as usual by making "minor" modifications to the vacuum 
solution. 

Another strange feature is the non-selfconjugacy of the solutions, at least in general. This 
is a consequence of the non-hermiticity of the Hamiltonian (I2.20[) .R Such solutions have been 
investigated already in the study |EFS05| of the gl(2|l) spin chain, but also encountered in 
different contexts dealing with non-hermitian Hamiltonians |SWKOO] . Instead, all lowest 
lying degree 1 solutions are invariant w.r.t. a modified conjugation symmetry 



gl(2m|l) : = {4"+^-'=)}^^"" (5.8) 

.l(2m + I|l) : {xf*}tl = {:E("+^-'-'}^lt"" , k^m + l. (5.9) 



The symmetry (|5.9p looks surprising, because the Dynkin diagram in fig. [5] is not symmetric 
w.r.t. the transposition ak ^ an+i-k- However, notice that the roots k ^ 

m + I together with the particle hole transformed roots {xj'"^^''}^!^^^^' also solve the BAE 
corresponding to the transposed Dynkin diagram of fig. El see |Tsu98] for details. For zy(™+i) 
even, the solutions are (exceptionally) selfconjugate 



V 



even : {xf*}fl\ = {xf , fc = I, . . . , n + I . (5.10) 



The following picture of the distribution of roots in the complex plane in the thermody- 
namic limit is emerging from our numerical analysis 

gl(2Hl): <^ ' ' ' ' 'i.^"'^^ (5-11 
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L -^3 ~ -^3 2 ' ^3 ~ 3 

The gl(2m|l) complex root configurations where called ±-strange strings in |EFS05j . 

We see that the lowest degree 1 excitation looks pretty complicated. What about other 
degree 1 excitations? We have verified that properties (|5.8H5.9p hold for all low lying exci- 
tations with the same root numbers as in fig. [6] on the right and, more then that, for all low 
lying solutions with symmetric configurations of the root numbers v'^^'^ ~ Many of 

these excitations seem to tend to the form (|5.111 15.121) in the thermodynamic limit. These 
are the best understood excitations for which we shall give a continuum description in the 
next section. We have also noticed that among the low lying excitations with even, 
there is always a subset of selfconjugate solutions of the type (j5.10p . These have interesting 
properties, as we shall see in a moment. 

To complete the overview of low lying degree 1 excitations, let us mention that we have 
also identified solutions with symmetrical configurations of the root numbers, which do not 



^See IVla86| for a proof of selfconjugacy in the case of hermitian Lie algebra spin chain Hamiltonians. 
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seem to tend to the form (|5.11[ I5.12p in the thermodynamic hmit. For these solutions 
we have nothing to say, although their study might prove crucial in constructing an S- 
matrix description of the continuum theory in large volume. The situation is even worse for 
solutions with non symmetrical configurations of the root numbers, because ±-strange string 
configurations cannot be clearly defined. 

5.2.3 Selfconjugate solutions 

As we have said, for even and i^^''^ = , there is always a subset of selfconjugate 

solutions of the type (jS.lOl) . It is straightforward to show that the gl(2m|l) BAE (|2.10p 
subject to the constraints (j5. 8115. 1(1)) are equivalent to the BAE of the so(2m + 1) fundamental 
spin chain of length L with periodic boundary conditions for the first m — 1 types of roots 
and antiperiodic boundary conditions for the spinorial roots. This is the generalization of the 
Takhtajan-Babujian subsector of the gl(2|l) spin chain of [EFS05| . The eigenvalues of the 
integrable Hamiltonian for this so(2m + 1) spin chain are given by eq. (|2.18p subject to the 
constraint (|5.10p . While this correspondence between solutions of gl(2m|l) and so(2m + 1) 
BAE is quite curious, the remarkable thing is that it extends to a correspondence between 
weakly excited solutions on both sides. In particular, the vacuum solution of the so(2to + 
1) fundamental chain in the thermodynamic limit |Mar91| (plus holes) is compatible with 
the form (|5.1ip when the constraint (IS.lOp is taken into account. We shall see how this 
correspondence can be used effectively in the next section. 

The situation is quite different for gl(2TO+ 1|1). Notice that assumptions Al— A2 and Rl 
of secHhold for solutions of gl(2m+l|l) BAE (P?TU|) satisfying ([010U| and, therefore, these 
can be restricted to solutions of gl(2m) BAE. However, the latter are very special, because 
in addition to satisfying the gl(2TO) BAE they must also admit multiple gl(2m + 1|1) lifts. 
According to the claim 2 of sec. [H gl(2m + 1|1) lifted Bethe vectors correspond exclusively 
to gl(2TO) solutions that do not admit multiple lifts. Presumably, what happens is that the 
additional conditions satisfied by these gl(2m) solutions obtained by restriction ensure the 
vanishing of the corresponding gl(2TO) Bethe vectors. We have checked that there are many 
other low lying solutions satisfying the constraints (|5.9[ 15.10]) . However, it is most confusing 
that many low lying solutions that do not satisfy the constraints (|5.91 15.10p in finite volume 
seem to have a thermodynamic limit (j5.12p that does satisfy the constraints (|5.9l 15.101) . We 
are not sure how to interpret this behavior, although we are tempted to believe this indicates 
that the large volume limit is a subtle issue. Therefore, we concentrate in the next section 
solely on the continuum description of gl(2m|l) degree 1 excitations of the type (|5.81 15.11]) . 

6 Continuous limit 

In this section we shall consider the spectrum of the spin chain of sec. [5] in the thermody- 
namic limit L ^ oo. The homogeneous spin chains are gapless and the only thermodynamic 
quantity, as far as the spectrum is concerned, is the vacuum energy per site. The latter 
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does not provide any insight into the the continuous hmit of the chain, which is expected to 
be governed by a conformal field theory (CFT). Probing this CFT requires computing scal- 
ing corrections to the spectrum jCar86j . However, these are much harder to compute, e.g. 
see [KP91| . then thermodynamic quantities. Therefore, gaining insight into the CFT solely 
from the lattice is quite non-trivial. To avoid such complications, one can introduce a smooth 
gap in the spin chain. The continuous limit is then expected to be a massive deformation 
of the CFT. Many interesting quantities, such as the /3-function, the particle spectrum, S- 
matrices, can now be computed in the thermodynamic limit. The standard way to generate 
a gap in the homogeneous spin chains is by introducing an alternating inhomogeneity A in 
the monodromies (j2.21 12. 3p . called staggering |Fad96) . The source terms of the homogeneous 
BAE (PIB then change to 



\xY-il2j Xx'-p - K/2-il2) \x'f> + A/2-1/2 J 
while the eigenvalues of the new energy and momentum operators become 



E ■ 



(6.1) 



-.(1) . . , . . , ,M) 



2 ^ \ ^ 2 \ ^ 2 2 ^ ^ \ J 2 \ ^ 2 

where a\ = (— l)'^', (Tm+n = (—1)1*^+^1 and we have assumed L to be even. These expres- 
sions reduce the the previous ones l2.lSl in the limit A ^ 0. 

In this section we define and study the continuous limit of the staggered gl(n -I- N\N) 
chains of sec. [H First we show that the bosonic lift of the staggered chain is described 
in the continuous hmit by the gl(n) Gross-Neveu (ON) model [AL80alllT8QbllAL79) . This 
strongly suggests that the continuous limit of the staggered spin chain coincides with the 
gl(n -f iV|A^) GN model, because the latter contains the gl(n) GN model as a cohoniological 
subsector |CCMS10] . The identification of the continuous limit is the main result of the 
paper. In the second subsection we explore the particle content of the gl(2TO|l) GN model 
which docs not lie in the gl(2m — 1) GN model lift, that is the bosonic lift. 

6.1 Bosonic lift 

As we have explained in sec. 2] and 15.11 the solutions of gl(n) BAE corresponding to non- 
vanishing Bethe vectors admit a unique lift to solutions of gl(rt -f A^jiV) BAE and Bethe 
vectors. We called this subsector of the gl(n -|- N\N) spin chain the bosonic lift. The 
uniqueness of the lift implies that the fermionic roots of the lifted solutions of gl(n -I- A^jA^) 
BAE, in the form schematically represented in fig. [SI are uniquely determined by the even 
roots. By the restriction procedure of sec. HI the latter solve gl(n) BAE. Therefore, all 
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1 • • ■ m 

0-<X^^ so(2m + l)l^®^ 

O 0# osp(1|2to-2)T/®-^ 

1 • • • m — 1 

Figure 7: Dynkin diagrams represent BAE for the corresponding chains. Double arrows 
indicate the equality of the corresponding set of roots. 

dynamical degrees of freedom (associated to holes) are determined by the even roots. The odd 
roots behave as auxiliary quantities useful for defining lifted Bethe vectors. So, we conclude 
that the BAE of the gl(n) spin chain describe entirely the bosonic lift of the gl(7i + N\N) 
spin chain. We recall some of the old results on the continuum limit of gl(n) spin chains and 
then reinterpret them in the gl(n + N\N) context. 

6.1.1 Vacuum energy 

According to sec. 15. 1[ the vacuum energy of the gl(n + N\N) spin chain coincides exactly 
with the vacuum energy of the g\{n) spin chain. The vacuum energy Coo = limL->.oo of 
the homogeneous gl(n) spin chain V^^ was computed in |Sut75] . This calculation is easily 
generalized for the staggered g\{n) spin chains {V Cg) V*)^^ 

+ > (6-2) 

where ip is the digamma function. Inserting n — 1 and A = we get the result eoo = —4 
of [EFSOSj . For n — 2ra — 1 the vacuum energy coincides exactly with the vacuum energy 
of the staggered so(2to + 1) and osp(ll2m — 2) fundamental chains of length L. These 
coincidences where probably first noticed in |SWK02llMar95j . The correspondence between 
the spectra of gl(2m|l), gl(2m — 1), so(2to + 1) and osp(l|2m — 2) chains is easily understood 
at the level of BAE from fig. [T] The non-trivial feature of the correspondence is that it relates 
low lying with low lying excitations. It is exactly this type of correspondences that lie at the 
heart of the dualities observed in [SPlOj . 

6.1.2 Particle spectrum 

Let us first recall some of the results on the continuous limit of gl(n) chains with real vac- 
uum solution. Going on the lines of jALSOb) . the gl(n) staggered spin chain {V ® V*)®'" 
can be shown to posses n — 1 gapped excitations branches associated with the holes in the 
distributions of n — 1 types of roots. For A big, the gap is of order e^^"^/". The continuous 
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limit is performed by introducing a physical lattice spacing a and then by taking the limit 
a 0, A — > oo such that 

m = a'^e-'''^/" (6.3) 

remains fixed, defining the mass scale of the theory. If we correct for the Fermi velocity 
c = At: /n, then in the continuous limit each excitation branch yields a relativistic dispersion 
relation with mass 

mfe = 2rnsin — , fc = l,...,n — 1. (6.4) 

n 

Therefore, we interpret it as a particle in some relativistic quantum field theory. The lowest 
lying excitation of the spin chain corresponding to a traceless tensor to^X, fi) of admissible 
shape (|5.ip is reproduced by a real solution of BAE with root and hole numbers fixed by the 
highest weight of the tensor |DN98j . This excitation is associated to a scattering eigenstate 
in the field theory of energy 

Ex.^, - -Bvac = ^ my^ cosh Oi + ^ m^; cosh 9i , (6.5) 

i=l 1=1 



where 6, 9 are the physical rapidities parametrizing the hole positions. Eq. (j6.5p suggests that 
the fc-th particle (j6.4p corresponds to either a covariant l*"' or contravariant P antisymmetric 
tensor. This is confirmed by an ^-matrix calculation |AL79[[DN98| . The end result is the 
S"- matrix of the gl(n) Gross- Neveu model [KKS79j . For instance, 1 and I particles scatter 
with themselves with the S'-matrix 

5n(0) = Snie) = Zn{9)^±^l^ , S,, = Z,m (l ^ „ ,^ ) , (6.6) 

1+1710/271 \ 11/2 + 1110/271 J 

where P and Q are permutation and Temperley-Lieb operators defined in sec. [2] and 

r(l-l/n-»0/27r)r(l + zg/2^) _ r(l/2 - ie/27T)T{l/2 - l/n + %e/2i:) 

^ r(l - l/n + i6i/27r)r(l -i6l/27r) ' ^ r(l/2 + z6'/27r)r(l/2 - l/n - i6'/27r) 

are, so called, dressing factors. The GN model coupling is related to the staggering g = A^^ 
and eq. (|6.3I) yields its /?-function. 

There is one important remark. Multiparticle states have zero total u(l) = gl("-)/ sl{n) 
charge, because the spin chain {V ®V*)'^'" has zero u(l) charge. Therefore, it is not possible 
to fix the u(l) charge of the individual particles in the spin chain. The allowable multiparticle 
configurations in the spin chain make the value of this charge irrelevant. However, the field 
theories with u(l) charged and uncharged particles (j6.4p are clearly different: the first one 
has 2n — 2 irreducible multiplets of particles, while the latter, which is the GN model, has 
71 — 1 irreducible multiplets. 

We now reinterpret the results of the above calculations in the framework of the bosonic 
lift of the gl(n+A^| A^) spin chain. The continuous limit of the gl(ri-|-A^| A^) spin chain is defined 
by the same eq. (j6.3|) . It already implies that the /3-function of the continuous theory does 
not depend on A^. The n — 1 excitation branches of the bosonic lift are interpreted again as 
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relativistic particles of mass (|6.4p . Eq. (|6.5p gives now the energy of a scattering eigenstate 
corresponding to a gl(n + N\N) traceless tensor to{X,fi) of g\{n) admissible shape (|5.1I) . 
However now, it suggests that the bosonic lift contains 2n — 2 distinct multiplets of particles 
corresponding to the gl(n + TVjA^) covariant antisymmetric tensors 1,1^,..., and their 
antiparticles, that is the contravariant antisymmetric tensors 1, P, . . . , 1"^^. Notice that the 
number of particles does not depend on whether they carry or not a u(l) charge, which again 
cannot be fixed in the spin chain. This is because the sl(n + |iV) antisymmetric tensors 1*^ 
and 1"^*^ arc not isomorphic anymore. 

The suggested gl{n+N\N) symmetry of the particles in the bosonic lift must be confirmed 
by an S'-matrix calculation. While we do not know the full ^-matrix, one can compute its 
restriction to the bosonic lift. The calculation is formally the same as for the gl(n) spin 
chain, because the BAE are the same. Thus, the eigenvalues of the gl(n + N\N) spin chain 
S'-matrix restricted to the bosonic lift coincide with the eigenvalues of the gl(7i) spin chain S'- 
matrix. More precisely, the eigenvalue of the gl(n-|-iV|7V) spin chain S-matrix on a scattering 
eigenstate corresponding to a gl(n -I- N\N) traceless tensor to (A, of gl(rt)-admissible shape 
does not depend on N. Due to the fact that the Q-cohomology of a tensor product of 
gl(n-f iV| A^) representations is the tensor product of their Q-cohomologies [CCMSIO] all these 
eigenvalues are compatible with our earlier assumption that the particles of the bosonic lift 
are co- and contravariant antisymmetric tensors of gl(n)-admissible shape. 

Finally, the scattering of particles in the bosonic lift can be computed exactly when no 
scattering eigenstates outside the bosonic lift are generated. Thus, particles 1 and 1 scatter 
with the same S-matrix (|6.6p . where this time P and Q act on tensor products of gl(n-f A^|A^) 
representations. Notice that we did not assume any crossing or unitarity to derive Sij from 
Sii. The question now is what field theory can reproduce these S-matrices? Given its 
formal similarity with the gl(ri) GN S-matrix, the most obvious candidate is the gl{n + 
N\N) GN model. There are deeper reasons to believe this. First of all, the Q-cohomology 
of the gl(n -|- 7V|A^) GN model is the gl(n) GN model. In particular, their /3-functions 
coincide |CCMS10j . Secondly, there is no doubt that the perturbative 1 /n calculations of the 
gl(n) GN S-matrix can be generalized. While we did not carry out an honest perturbative 
S-matrix calculation as in |KKS79) , the result should certainly be (|6.6p . This is because 
the Feynman rules are g\{n + N\N) invariant tensors and their algebra, generated by tensor 
multiplication and contraction, is a representation of the walled Brauer algebra Bl^lIu) of 
sec. 13.11 and 13.21 Thus, as long as one computes an invariant tensor, such as the S-matrix, 
starting from other invariant tensors, such as the Feynman rules, the result cannot depend 
on N. Thirdly, the detailed analysis of |EFS05[rSS07j suggests that the continuum limit of 
the homogeneous gl(2|l) spin chain is the si(2|l)i WZNW model, that is the gl(2|l) GN 
model at zero coupling. Finally, we mention the similarity with the situation for osp(i?|2S) 
GN models |BL00i[SWK02llSP10| . Therefore, we conjecture that the continuous limit of the 
gl(n + N\N) spin chain {V ® F*)®-^ is the gl(n + N\N) GN model. 
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6.2 Degree one excitation 

In this section we investigate degree 1 excitations of the gl(2m|l) spin chain which are of 
the type (jS.Sp and tend to the form (jS.llD in the thermodynamic hmit. As explained in 
sec. 15.21 the lowest lying excitation of degree 1 is precisely of this type. It corresponds to the 
antisymmetric tensor to{l™, 1™) and its root numbers are represented in fig. El 
In the thermodynamic limit the roots 

{-f - ^..™..-) k = l,...,m-l (6.7) 

and the centers {a:^^''}^]^, of the ±-strange strings become dense on the real line. We now 
derive the Lieb equations satisfied by their densities. Collecting all BAE, in the grading 
[+™ - +"], we get 



-1 = eiix^ A/2)^/2g^(^ ^ K/2f'^E^^\x)E^^l{x) , x £ {xf^})t\ (6.8) 

-1 = i^iV'^ {^)E^2^ (x)Ei\+^^ (x) , fc = 2, . . . , TO - 2 , xe {xf'> }f^l 

2 2 2 2 J J 

+1 = E^T'\x)E^:^l~'^ ix)Ei^'> ix)E[^^ [xf , X e {xf ^}fii 

2 2 J J 

where e[^\x) = 11^=1 ^t{x — x'j'^), for k — 1, . . . ,to — 1,±. The first two sets of BAE 
equations is just a rewriting of eqs. (j2.10p . To get the last three we have multiplied the 
equations for same root appearing symmetrically in the Dynkin diagram |5] at positions m — 1 
and m + 2 or TO and to + 1. Notice the + sign appearing in the last two sets of equations. 

Define the densities of roots Pk{xf^) :— limi_>oo l/i(a;j^i_\ — Xj'^^) and holes p'^{x) := 
J2]jLi ^(x — Cj*^^ )/i for fc = 1, . . . , TO — 1, ±. Taking the logarithm of eqs. (|6.8p and derivating 
w.r.t. the spectral parameters we get the Lieb equations 

ioi * {Sa + S_a) = p'l + Pi + a2 * Pi — ai * P2 (6.9) 
= Pk + Pk + 0,2 * Pk - ai * (Pfe-i + Pfc+i) , fc = 2, . . . , TO - 2 

= Pm-l + Prn-l + a2 * p,„_l - fli * p,n~2 ~ ^{as + ai) * {p+ + p^) 

Q = p± + p± + a2 * p± + ^ai * p^ - (a 3 + ai ) * pm-i , 

where 5x is the Dirac distribution centered at x, at{x) — ^^^^j— = ^ ^i^/^i/^ and * denotes 
the convolution {f *g){x) — dy f{x ^ y)g{y). From eq. (j6.ip . the energy and momentum 
of a root configuration (|6.9p is 



E ^ -2TrL / dxpi{x)[ai{x + A/2) + ai{x-A/2)] 

J — OO 

P = L dxpi{x)[ei{x + A/2) + 9i{x-A/2)] 
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To solve for the root densities in terms of the hole densities it is useful to define the 
distributions 

_ p+ +p- _ p+ - 

Ps- ^ , Pa- ^ 

2 ' 2 ' 

Notice that the distributions pjj, pa split from the rest of eqs. (j6.9p 

= pjj + Pa + 02 * Pa - 2ai * Pa , (6.10) 

while the Lieb equations for the remaining densities pi, . . . ,pm-i,Ps look exactly like the 
Lieb equations for the staggered so(2to + 1) fundamental chain of length i, e.g. |Mar91) . The 
resolvent for the latter can be found in |ORW87l . Computing the solution of eqs. (|6.8p with 
no holes one gets 



2 sin ^ cosh TT 2 sin ^ cosh tt 



2x-K 



0£j Oi-i-i. VjV/Oi-A l\ ^ Cli-i-i. VjWOJ.i. IK , , 

Pk^^) n cosh HIL(H£±A) _ cos ^ ^ cosh HllHi^ - cos 2^ ' ^ ' 

n n n n 

''^"^ ^ ^^^^^[^^ ^ cosh ' ^ = l'---''"-l' 

n n 

where n — 2m — 1. We want to stress that this is a formal solution that does not correspond to 
any state in the spin chain. This is because the degree 1 excitations we are considering always 
contain holes. However, it is important to realize that the energy of the formal solution (|6.1ip 
coincides with the actual vacuum energy (|6.2p . We have already explained in sec. 16.1.1) why 
this happens. As usual, the energy and momenta of holes can be expressed in terms of hole 
less densities 

e,{0=^npl{0, ^^2npl{0 

^±(e)=2.p°(0, ^ = 

Restricting to low lying states |^| ^ A one gets the expected relativistic dispersion relations 

ek{d) — 2camkCOsh9 , Pk{d) = 2amkSmh9 (6.12) 

e±{9) — cam cosh 6* , P±{0) = am sinhfl (6.13) 

with the gl(fT.) Fermi velocity c = Att/ti and physical rapidity 6 = 27r^/n. The masses are the 
same as in eqs. (|6.3l 16. 4p . 

It is certainly reassuring that degree 1 excitations have the same mass scale m as degree 
excitations. The doubling of masses (|6.12p w.r.t. the masses of particles in the bosonic 
lift can be understood as follows. We have observed in sec. 15.2.21 that the imaginary parts 
of the odd roots tend to vanish for large L and, therefore, we have neglected them in the 
thermodynamic limit. However, due to the conjugation symmetry (15. 8p of the solutions, this 
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approximation leads to an identification (I6.7P of distinct roots. Tfierefore, tlie positions of 
holes in the distribution of even roots {x^j'''}'jLi and {a;^-^™^^ '^''ij'li ^ ' have also been 
identified. So, the excitation (j6.12p is in fact a state of two particles 1*^ and P each of mass 
nik and rapidity 9. These are the old particles from the bosonic lift. On the other hand, the 
particles (|6.13p are new. 

The artificial binding of holes (j6.12p induced by the thermodynamic limit does not allow 
to derive meaningful S'-matrices from eqs. (|6.9p Most probably, to correct the approach and 
recover the degrees of freedom lost in the thermodynamic limit one needs to distinguish 

(k) (k) 

between the two types of roots Ima;^- > and Imx^ < 0, exactly as we did for the strange 
ib-strings (jS.llI) . However, a finite volume treatment will be required, because the imaginary 
parts of even roots vanish in the thermodynamic limit. 

The need for a finite volume approach can also be seen from eq. (|6.10p . Solving for the 
Fourier transform of Pa{p) = J dp exp{—2Tripx) pa{x) one gets 

MP I ^^^\&y) = AeHpi-^^p(«'^'+«'-')5i!iM!^^iL!) , (6.14) 

FaKFl^ ) sinh^^ ^ ' 

where we have considered only a single pair of ±-holes. This is no loss of generality, because 
from (|6.9p the numbers of i-holes is always equal 

n± = 2(:/("-i) -N+- iV_) = 2(z^("-i' - j/(")) . 

and, therefore, they always come in pairs. The Fourier transform (j6.14p is singular and clearly 
must be regularized, because one has to satisfy the constraint 

h := lim = L / dxpa{x) = pa{0) , (6.15) 

where 6 is a fixed real number parametrizing the state. On a lattice of length L the the 
momentum can take a minimal value of p^^l/L. Using this value as a regulator one gets for 
a pair of ±-holes from eqs. (|6. 14116.15]) 

^ ^ 2i ' 

which clearly shows that for 6 ^ the deviations of strange i-strings from the form (|5.11l) of 
the solution in the thermodynamic limit have to be taken into account. Notice that the energy 
does not depend on the continuous parameter (j6.15p parametrizing the state. Therefore, it 
is tempting to conclude that there is a continuum of new particles (|6.13p . 



7 Conclusions and Outlook 

We have put on firm grounds the relationship between gl(n + N\N) integrable spin chains 
with n fixed. This allowed us to prove that all gl(n + N\N) spin chains {V ]/*)®^ with 
n, N > possess in the continuum limit 2n — 2 multiplets of massive particles which scatter 
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with gl(n) Gross-Neveu like 5'-matrices, namely their eigenvalues do not depend on N. We 
concluded that the continuum theory is the gl(M|A^) Gross-Neveu model. Evidence that 
the massive spectrum is much richer, possibly continuous, was established on the example 
of gl(2m|l) chain. Finally, our analysis of the thermodynamic limit strongly suggests that 
understanding the nature of new particles requires a finite volume treatment. 

The question that begs the quickest answer is how to close the fusion of S'-matrices (|6.6p of 
the gl(M|iV) Gross-Neveu model starting with just the vector multiplet and its antiparticles. 

The gl(Af|A^) spin chains require a separate treatment, which we hope to report on later. 
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A Highest weight vectors 

Let T be the indecomposable direct summand of t{X,fi) with traceless submodule to(A,/x). 
Without loss of generality, one can assume that (A, /i) is of admissible shape in the sense of 
eq. 13.41 The filtration (13. ip implies a similar filtration for T 

to(A,A*) =To cTi C---cTf=T 

where we have defined the submodules Tn — T O t„(A, /i). Again as for t(A, ^) one has 

TjTr,.,^ io(A',M'), (A.l) 

where the elements (A',/i') of D„ must be pairs of diagrams with f — n boxes such that 
A' C A and fi' C fi. In particular, they are also admissible. According to assumption 13.21 
every direct summand Iq^X',^') in eq. (IA.1|) can be generated from a highest weight vector 
v(A', /x')/r„_i e to{X',fi') w.r.t. some Borel subalgebra bs, where v{X',fj.') G r„. Assume 
that v(A', fi') is a highest weight vector w.r.t. bs. Then v(A', fi') cannot generate more then 
a Kac module within T„, which according to assumption 13.21 is precisely to^X',^'). Thus, 
to(A',/i') is a submodule of T„ and therefore also of T. This completes the proof of claim [T] 
Notice that if T(A, /i) is of the usual "diamond" form, that is its socle and top are irre- 
ducible, then the existence of two submodules to{X, iJ,),to{X' , fi') C T{X,iJ,) contradicts the 
simplicity of the socle and, therefore, all highest weight vectors of T(A, fi) belong to to{X, fi). 

B Cohomological reduction of Kac modules 

In this section we prove the claim [2] of sec. U and adopt all the notations leading to it. 
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Assume that we have a solution for p < P and a corresponding Bethe vector uti. It 
comes together with a solution for p = P and a corresponding Bethe vector uj2 ■ The transfer 
matrices (|2.4p yield the same eigenvalues on uji and u}2- Therefore, they must be part of a 
reducible indecomposable module or, more precisely due to corollary 1 of sec. 13. 2[ part of a 
Kac submodule of the latter. Moreover, because wt(a;2) < wt(aji) the module generated by 
UJ2 must be a submodule of the module generated by uji 

ll{g\{M\N)) ■ L02 C U{g\{M\N)) ■ Loi C /^,(A). (B.l) 

Here A is the highest weight of the Kac module w.r.t. some Borel subalgebra b. Obviously, 
the highest weight must be atypical (A,ckfc) = in order for (|B.ip to hold. Recall that 
the Borel subalgebras bs and b of gl(M|iV) define each a Z-graded decomposition of the 
latter |FSS96j 

gl(Af|iV)~g-i©0O ©fli , fli Cb 
gl(Af|7V) ^0s'©0l©0s , flsCbs. 
Recall the definition of Kac modules 

K,{k) ^ /\ (g-i) ®z,(gOeBi) 5o(A) , (B.2) 

where Sq (A) is an irreducible representation of go of highest weight A trivially extended to a 
representation of 0° © g^. Then, due to the embeddings (jB.ip . one must have Q G 0^^ n flj^^- 
Notice now that Kac modules have zero superdimension. This is a necessary condition for 
vanishing Q-cohomology. Indeed, given Q G 0^^ it is straightforward to show that the Q- 
cohomology of the Kac module defined as in (IB.2p vanishes. For more details, see the proof 
of [CCMSIO] for the vanishing of the cohomology of projective modules. 

Notice that while projective modules have zero cohomology w.r.t. any choice of Q, the 
Kac module (|B.2p . on the other hand, has zero cohomology only w.r.t. Q £ 0^^. 
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